Benard-von Karman vortex street in an exciton-polariton superfluid 
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The dynamics of an exciton-polariton superfluid resonantly injected into a semiconductor micro- 
cavity are investigated numerically. The results reveal that a Benard-von Karman vortex street is 
generated in the wake behind an obstacle potential, in addition to the generation of quantized vortex 
dipoles and dark solitons. The vortex street is shown to be robust against a disorder potential in a 
sample and it can be observed even in time-integrated measurements. 
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I. INTRODUCTION 

A cylindrical obstacle moving in a classical viscous 
fluid generates a vortex- antivortex pair in its wake. At 
large Reynolds numbers, such vortices develop into a pe- 
riodic pattern known as the Benard-von Karman vortex 
(BvK) street^ Such dynamics are also found in super- 
fluids for which vortex circulation is quantized and vis- 
cosity is absent. It was numerically demonstrated 3 that 
quantized vortex- antivortex pairs, which we call vortex 
dipoles, are shed in the wake of an obstacle moving in a 
planar superfluid above a critical velocity. This theoreti- 
cal prediction has been realized in Bose-Einstein conden- 
sates (BECs) of atomic gases, in which quantized vortex 
dipoles have been created using an moving obstacle po- 
tential produced by a laser beam^£ BvK vortex streets 
in superfluids (i.e., periodic and alternating creation of 
quantized vortices and antivortices that form a vortex 
street with a long lifetime) have also been predicted for 
atomic BECs^ but they have yet to be demonstrated 
experimentally. 

Exciton-polariton superfluid flow past an obstacle has 
recently been demonstrated using a semiconductor mi- 
crocavity^ - — In these experiments, polaritons with a 
controlled in-plane momentum are coherently injected 
into a sample by a pumping laser, and an obstacle po- 
tential is produced by a defect in the microcavity^" 
by a continuous- wave (cw) laser field, 11 or by etching 
the sample^ 2 A C her enkov- like pattern and oblique dark 
solitons are observed for supersonic flow of polariton con- 
densates^ Oblique dark solitons formed behind the ob- 
stacle decay into quantized vortices For subsonic flow, 
quantized vortex dipoles are produced in the wake^~— 
These phenomena have been studied theoretically by sev- 
eral researchers ji 3 -"— 

A polariton condensate differs from superfluid helium 
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and an atomic gas BEC in that it is a nonequilibrium 
open system. The polaritons have a short lifetime of ~ 
10 ps, which is comparable to the time scale of the rele- 
vant dynamics. In the experiments described in Refs.[8||9|, 
polaritons are constantly replenished by pumping with a 
cw laser and the system reaches a nonequilibrium steady 
state, in which the pumping balances the loss. Another 
difference from atomic systems is that the polariton is a 
coherent superposition of a quantum-well exciton and a 
cavity photon and the interparticle interaction originates 
only from the former; in other words, a polariton super- 
fluid is a nontrivial two-component system. Because of 
these differences, it is by no means obvious that a polari- 
ton superfluid shares the same dynamic phenomena as 
superfluid helium and atomic gas BEC. 

In this paper, we investigate the dynamics of an 
exciton-polariton superfluid passing an obstacle potential 
and show that a superfluid BvK vortex street 6 emerges in 
this system, as well as vortex dipoles and dark solitons. 
We also show that the vortex street is not destroyed by 
a disorder potential, which is present in a realistic sam- 
ple. We propose a measurement method to identify a 
vortex street by time-integrated imaging for the present 
cw pumped system. 

This paper is organized as follows. Section [II] formu- 
lates the problem, Sec. IIIII presents the numerical results, 
and Sec. HVl gives the conclusions of the study. 



II. FORMULATION OF THE PROBLEM 

We consider a system of quantum-well excitons and 
cavity photons in the mean-field theory. The mean-field 
wave functions of excitons ipx and photons ipc are as- 
sumed to obey the two-component nonlinear Schrodinger 
equation in two dimensions, 13 
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The polariton Hamiltonian in Eq. (pQ) is given by 
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where the diagonal elements are the dispersion relations 
of an exciton and a photon, and £7r is the Rabi fre- 
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quency of the exciton-photon coupling. In Eq. (pQ), 9 is 
the exciton-exciton interaction coefficient, 7x and 7c are 
respectively the exciton and photon decay rates, and V 
is the potential for cavity photons. The last term on the 
right-hand side of Eq. (pQ) describes the coherent quasi- 
resonant pumping of photons by an external laser beam 
with a spatial profile F, an in-plane wave vector fe p , and 
a frequency uj p . 

Diagonalizing the noninteracting Hamiltonian in 
Eq. ([2]), we obtain the eigenfrequencies of the upper and 
lower free polaritons as 

ux + cj c ± VVx -u c ) 2 +4^| 
^± 2 ' ^ ' 

The pumping frequency uj p is close to cj_ to resonantly 
excite the lower polaritons. We define the detuning as 
6 = uj p — U-(k p ). The exciton is assumed to have the flat 
dispersion, 0Jx(k) = c<; x , and the dispersion of the cavity 
photon is approximated to be uc(k) = uoq + ftk 2 / (2mc), 
where mc is the effective mass of a cavity photon. In the 
following, we assume u x = Uq. 

We numerically solve Eq. (pp) using the pseudospectral 
method, 1 - which imposes periodic boundary conditions. 
The spatial range in the numerical calculation is taken 
to be sufficiently large that the boundary conditions do 
not affect the results. The initial state is the vacuum 
state of excitons and photons, = = 0, plus small 
white noise to break the numerically exact symmetry. 
The results do not qualitatively depend on the detail of 
the noise. 



III. NUMERICAL RESULTS 

In the following numerical calculations, we use mc = 
2x 10 _5 ra e , where m e is the electron mass, hflR = 5 meV, 
and g = 0.01 meV/im 2 , and we assume 7x = 7c = 7- 
The pumping function F(r) is assumed to have the form 

J?M - I F ° ( ax < 1^1 < y v and x > °)' (A\ 
r{r) -\0 (otherwise), [V 

which excites the polaritons in the region of the two tri- 
angles shown in Fig.[T](a). We restrict the pumping area 
as in Eq. ((4]), since if the whole space is pumped, the 
phase will be locked and no vortices will be generated. 
Similar pumping schemes have been employed in theo- 
retical^ and experiment alii studies. The parameters in 
Eq. (|4]) are taken to be HFq = 38.2 meV, a = 0.4, and 
y v = 50 fim. The in-plane wave vector of the pumping 
beam is fe p = k p e Xl where e x is the unit vector in the 
x direction. Polaritons pumped in the triangular areas 
thus have momentum in the x direction and flow into the 
region between the triangles, as shown in Fig. [1] (a). An 
obstacle potential is assumed to be a Gaussian potential 
as 




FIG. 1: (Color online) Density l^cl 2 and phase arg(^c) pro- 
files of the photon wave function at t = 400 ps. (a) k p = 0.25 
/im" 1 without an obstacle potential, (b) k p = 0.25 /im" 1 and 
d — 2 /im, (c) k p — 0.25 /xm _1 and d = 3 /im, (d) k p — 0.25 
/im -1 and d — 5 /im, and (e) k p — 0.5 /im -1 and d — 3 /im. 
The solid circles in (b) and (c) respectively indicate vortex 
dipoles and co-rotating twin vortices, which are released from 
the obstacle potential. The dashed circles indicate the loca- 
tions of the obstacle potential and the filled blue circles and 
red squares indicate clockwise and counterclockwise vortices, 
respectively. The detuning is Ti5 — 0.7 meV and the lifetime 
is 7 _1 = 30 ps. The field of view of each panel is 150 x 110 /im 
and the origin is located at the center of the left edge of the 
panel. See Supplemental Material for movies of the dynamics 
in (b)-(e). 



V(r) = V exp [-(* - x ) 2 /d 2 - y 2 /d 2 } (5) 
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with Vo = 38.2 meV and xo = 40 /im, which is located be- 
tween the triangles. Thus, polar itons flowing between the 
triangular areas hit the obstacle potential and generate a 
wake in the x direction. The healing length and the Bo- 
goliubov sound speed near the obstacle potential are re- 
spectively estimated to be £ ~ ^/(mLP^LP^Lp) 1 ^ 2 — 1.5 
/im and v s c± (^lp^lp/^lp) 1//2 — 1-8 x 10 6 m/s, where 
^lp, #lp, and tilp are respectively the mass, the effec- 
tive coupling constant^ and the density of the lower 
polariton. 

Figure Q] shows profiles of the photon wave function 
ipc at t = 400 ps. In Fig. [T](b), quantized vortex dipoles 
are generated behind the obstacle potential. After the 
vortex dipoles are released from the potential, they al- 
ternately align in the wake, as shown in Fig. Q](b). Such 
alternate alignment of vortex dipoles in superfluids is also 
observed in Refs. 6.18, Figure [H (c) shows the situation 
for a larger obstacle potential. In this case, two vortices 
with the same circulation form a pair and clockwise and 
counterclockwise pairs are shed alternately. This vor- 
tex shedding dynamics is very similar to that in Ref. 0, 
which is a superfluid analogue of the BvK vortex street. 
This periodic vortex configuration has a very long life- 
time 6 when there is no dissipation or losses. Thus, a 
small obstacle generates vortex dipoles and a large ob- 
stacle generates a BvK vortex street, which is consistent 
with the case of an atomic BEC (Fig. 3 of Ref. 6|). For a 
larger obstacle potential, vortex shedding becomes irreg- 
ular as shown in Fig. [T] (d). The flow velocity in Figs. [T] 
(a)-[T] (d) is Tikp/rriLP ~ 0.4v s , which is comparable to 
the critical velocity for vortex nucleation obtained from 
the Gross-Pitaevskii equation. 3 Figure Q] (e) shows the 
case for faster flow with flow velocity of ~ 0.7v s . Dark 
solitons and vortex dipoles are generated in the wake, as 
experimentally observed^ The vortex patterns shown 
in Fig. [I] can be observed for a shorter polariton lifetime 
(e.g., 15 ps) if the polariton density near the obstacle is 
kept large by increasing the pump intensity or narrowing 
the space between the triangular pumped regions (data 
not shown). For all the cases in Fig. [TJ the exciton wave 
function has a similar profile to that of the photon 
wave function ipc- 

In microcavity samples used in experiments, a disorder 
potential is inevitable; it is of the order of 0.1 meV, even 
for a good sample. 19 Figure [2] shows the effect of the 
disorder potential on vortex street formation, where the 
potential shown in Fig. [2](b) is added to V(r) in Eq. (pQ) 
while the other parameters remain the same as those in 
Fig. [11(c). Figure [2] (a) shows that the BvK vortex street 
is robust against a disorder potential in a realistic sample. 

In the present situation, it is impossible to per- 
form time-resolved measurements, as were performed in 
Refs. l8l l 1 Oil 12l For time-resolved measurements, polar itons 
must be pumped by a pulsed laser beam, which is split 
and used as a reference in the interferometer with a de- 
lay time, enabling time-resolved images to be obtained by 
performing repeated measurements. On the other hand, 
polaritons are pumped by a cw laser in the present system 



(a) (b) 
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FIG. 2: (Color online) (a) Density profile l^cl 2 in the pres- 
ence of the disorder potential shown in (b). The disorder po- 
tential is generated by setting random numbers on each site 
and cutting off short- wavelength Fourier components. The 
parameters are the same as those in Fig.[T](c). The black and 
white circles indicate clockwise and counterclockwise vortices, 
respectively. See the Supplemental Material for a movie of the 
dynamics. 
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FIG. 3: (Color online) (a) Time-integrated density of the 
photon wave function J(r) = f \ipc(r, t)\ 2 dt for the dynamics 
in Fig.[T](b). (b) I for the dynamics in Fig.[T](c). (c) Density 
\G\ 2 (left) and phase arg(C) (right) of the first-order coher- 
ence G(r) — J ipc( r vefi t)tpc(r, t)dt for the dynamics in Fig.[T] 
(b). (d) \G\ (left) and arg(G) (right) for the dynamics in 
Fig. [2(c). The time integrations in / and G are taken between 
t = 400 and 1400 ps (the patterns are independent of the 
upper integration limit if the time integration is sufficiently 
long). The position of the reference is (x re f,2/ref) = (46,4), 
which is marked by the crosses. The field of view of each 
panel is 150 x 110 /im and the origin is located at the center 
of the left edge of the panel. 
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and hence only time-integrated images can be obtained. 
Figures [3] (a) and[3](b) show time-integrated density pro- 
files J \ifjc(r^t)\ 2 dt of the dynamics in Figs. Q] (b) and [I] 
(c), respectively. Although the traces of vortex flow are 
visible, individual vortices are smeared out and cannot 
be discerned*- 9 To overcome this problem, we examine 
the time-integrated spatial coherence given by 

G(r) = J r c (rreut)Mr,t)dt, (6) 

where r re f is the position of the reference light source. 
The reference light from a small area at r re f is enlarged 
and interference with the whole image is measured using 
an interferometer with a variable arm length. 9 

Figures [3] (c) and [3] (d) show the density and phase 
profiles of G for the dynamics in Figs. [T] (b) and [T] (c), 
respectively. The reference r re f is positioned obliquely 
behind the obstacle (crosses in Fig. [3]), where the den- 
sity oscillates with time due to periodic vortex shedding. 
Therefore, time integration of G is performed stroboscop- 
ically at the same frequency as the change in the vortex 
pattern, preserving the information of vortices, as shown 
in Figs. [3] (c) and El (d). Figure [3] (c), which corresponds 
to the vortex dipole generation in Fig. Q] (b), exhibits 
vortex dipoles in the phase profiles (left panel of Fig. [3] 
(c)). The phase profile in Fig.[3](d) also reflects the BvK 
vortex street in Fig. (U(c). Thus, the time-integrated co- 
herence G contains information on the vortex patterns 
and it enables us to identify vortex street formation^ 



IV. CONCLUSIONS 

We have investigated the dynamics of exciton- 
polariton superfluid flow passing an obstacle potential 
in a semiconductor microcavity. Numerically solving the 
two-component nonlinear Schrodinger equation in Eq. (pp) 
reveals a superfluid BvK vortex street (Fig.[T](c)), as well 
as alternately aligned vortex dipoles (Fig.[T](b)) and dark 
solitons (Fig. Q] (d)). The formation of the vortex street 
is robust against a disorder potential, which is naturally 
present in samples (Fig. [2]). These periodic vortex pat- 
terns are time-dependent steady states that are attained 
by cw pumping of polaritons and time-resolved imaging 
cannot be performed. We showed that the spatial co- 
herence in Eq. (j6j) with an appropriate reference point 
r re f reflects the vortex patterns even though it is time 
integrated (Fig. [3]). We have thus demonstrated that a 
superfluid BvK vortex street can be generated and de- 
tected in an exciton-polariton condensate using current 
experimental techniques. 
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